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Abstract 

The (A/ = 1) decay D*+ D+ir^ and (A/ = 0) decay D*^ D+'j have been reconstructed 
using 90.4 fb^^ of data recorded by the BABAR detector at the PEP-II asymmetric-energy e'^e~ 
collider. The relative branching fraction of isospin violating decay to isospin conserving decay 
has been measured. Our preliminary result for this ratio is T{D*'^ Dfn^)/T{D*'^ —>■ Df^) = 
0.0621 ± 0.0049 (stat) ± 0.0063 (syst). In addition we reconstruct the decays D*^ D°it° and 
D*^ — > L)^7 for completeness and measure a preliminary relative branching fraction to be T{D*^ — > 
D'^Tr°)/T{D*^ L»°7) = 1.740±0.020(stat)±0.125(syst). Both measurements represent signihcant 
improvements over present world averages. 
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The hadronic decay^ D*^ — > Dfn^ violates isospin conservation and should be suppressed with 
respect to the radiative decay — s- -0^7- Cho and Wise ^1^ have calculated the decay rate for 
D*'^ — > -D^tt" using chiral perturbation theory; they describe the process as an isospin conserving 
decay D*~^ — > Dfr] followed by the conversion of the virtual rj meson into a vr'' meson through 
mixing. The radiative process D*~^ — > Df^ proceeds via magnetic transitions to which both the 
coupling to the charm quark and the light quark contribute. Since a neutral pion or kaon do 
not couple to the photon, only virtual states with a charged pion or kaon contribute to one-loop 
corrections of the magnetic moments. The difference in the kaon and pion masses naturally leads 
to SU(3) flavor-symmetry breaking. 

A previous measurement of the ratio T{D*+ L>+7r°)/r(L>*+ Df-/) by the CLEO Collabo- 
ration obtained a result of 0.062 + [j;[j?^(stat) ± 0.022(syst) using 14.7 t D+tt^ events [2]. 

The analysis we present here, using 90.4 fb~^ of data recorded by the BABAR detector at the 
PEP-II asymmetric e~^e~ storage ring, represents a significant improvement in precision and is 
particularly relevant given the recent observations of two new meson states [311]. The data 
have been recorded at and approximately 40MeV below the T{4S) resonance. 

The BABAR detector is described in detail elsewhere Charged particles are detected and 
their momenta measured by a silicon vertex tracker (SVT) consisting of five layers of double-sided 
silicon strip sensors and a cylindrical 40-layer drift chamber (DCH), both operating within a 1.5 T 
solenoidal magnetic field. Charged particle identification is provided by energy loss measurements 
in the SVT and DCH and by light recorded in an internally reflecting ring imaging Cherenkov 
detector (DIRC). Neutral particles are identified and their energies measured by an electromagnetic 
calorimeter (EMC) composed of 6580 CsI(Tl) crystals. 

mesons are reconstructed using the decay — > (f)7r~^ , (p ~^ K^K^ . Kaons are identified 
by combining the energy deposited in the SVT and DCH with the DIRC information. Tracks not 
identified as kaons according to the PID criteria are considered as pions. 

Our selection criteria are chosen such that the signal significance S'^/{S + B) of the channel 
D*'^ D'^TT^ is maximized. S and B refer to the expected numbers of signal and background 
events, respectively, which are obtained from a sample of Monte Carlo events. 

All combinations of K~^K~tt~^ candidates have been required to successfully fit to a common 
vertex; only combinations with a K^K'^ mass within 8MeV/c^ from the nominal (j) mass pj have 
been retained. 

The scaled momentum is defined as Xp{Df) = p* (Df) / p^^^{D^) , where p*{Df) is the momen- 
tum of the candidates in the center-of-mass frame, and p*na.x{Df) = \J E^^^^ — m{Dt)^ is its 
maximum value. Combinations of (j) and vr"*" candidates are retained if the scaled momentum is 
0.6 or greater. Since the kinematics of the cc fragmentation process produce charmed mesons with 
high momenta, this reduces the combinatorial background. 

We exploit the longitudinal polarization of the (j) meson to reduce continuum background by 
requiring that the absolute value of the cosine of the helicity angle 9h , which is defined as the angle 
between the (j) momentum direction in the rest frame and the momentum direction of one of 
the kaons in the (j) rest frame, is 0.3 or greater. 

The K^K^TT^ invariant mass distribution is shown in Fig. ^ By fitting the sum of a double 
Gaussian function (to represent the signal) and a third-order polynomial (to represent the back- 
ground) to this distribution, we find 73 500 it 300 events. Df candidates are retained if their 
invariant mass is within 12 MeV/c^ fr om the nominal Df mass jH] • 

^Inclusion of the charge conjugated state is assumed throughout this paper. 
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Figure 1: K^K^ir^ mass distribution. The dots represent data points. The sohd line shows the 
fitted function. The dashed hue indicates the portion of the fit associated with the background, 
candidates are selected from the region between the vertical dotted lines. 



Each vr candidate is reconstructed by combining two photon candidates. Each photon candidate 
consists of a calorimeter cluster which is not associated with a charged track and has an energy in 
the laboratory frame of at least 45 MeV. Additionally, to help remove the background from hadrons, 
we require the lateral moment jTj , which describes the shape of the electromagnetic shower in the 
calorimeter, to be less than 0.55. 

The TT^ candidates are retained if they have a momentum p* in the center-of-mass frame greater 
than 150MeV/c. Furthermore, the absolute value of the cosine of the decay angle 6*{tt^), which 
is defined as the angle between the direction of one of the photons in the vr'' rest frame and the 
direction of the tt" candidate in the center-of-mass frame, is required to be less than 0.85. For 
the isotropic tt" decay, the | cos6'*(7r'')| distribution is flat, while it peaks near 1 for random 77 
combinations. 

Only 77 pairs inside a mass window are retained. This window is defined as the region where 
a function fitted to the 77 mass distribution of true vr'^ — > 77 Monte Carlo events exceeds 0.2 its 
maximum value, which accommodates the asymmetric shape of the 77 mass distribution and takes 
different detector calibrations into account. A kinematic fit is applied to the surviving 77 pairs, 
constrained to the nominal vr'^ mass. 

After combining the Df and vr'^ candidates to reconstruct the decay Dfn^, we have 

fitted a function to the mass difference Am{Df7r^) = m(K^ K^tt^tt^) — m{K~^ K^tt^). The 
function is the sum of a double Gaussian function to represent the signal and 

/i(Am) = iV • - exp ^ . ^Am^ + a Am + b) (1) 

(where m[7:^) is the vr'^ mass, while A^, fi, a, and b are free fit parameters) to represent the 
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Am [GeV/c^] Am [GeV/c^] 

Figure 2: D*+ signals: (a) m{K+ K-ir+ir^) - m{K+ K-7r+); (b) m{K+ R-ir+'y) - m{K+ R- 11+) . 
The dots represent data points. The sohd hne shows the fitted function. The dashed hne indicates 
the portion of the fit associated with the background. 



background. The exponential term in the background function models the kinematic threshold; it 
is close to 1 in the signal region. 

The result of this fit is shown in Fig. I21^a). We find a signal event yield of 560 it 40. 

For the reconstruction of the decay D*'^ — > Dfj, calorimeter clusters that are not associated 
with a charged track are considered photon candidates if they fulfill the following requirements. 
The energy must be 50MeV or greater in the laboratory frame and 100 MeV or greater in the 
center-of-mass frame. The lateral moment must be less than 0.8. To reduce the background of 
photons coming from vr*^ decays, a photon candidate is discarded if it forms a vr^ candidate with any 
other photon candidate in the same event. Here, a 77 combination is considered a candidate 
if the invariant mass is in the range 115 < m(77) < 155 MeV/ and has an energy of at least 
200 MeV in the center-of-mass frame. 

To obtain the D*~^ — > Dfj signal event yield, we have fitted a function to the mass difference 
Am{D+j) =m{K+K-TT+-/)-m{K+K-Tr+). This function is a sum of a "Crystal Ball" function [Sj 




MAm)=N.{ \ J J ^/ - (2) 

- ) if (Am - ^)/cr > a 



(where A^, fi, a, n, and a are free fit parameters, and A and B are chosen such that the function and 
its first derivative are continuous at (Am — fi)/a = a) for the signal, and a third-order polynomial 
for the background. The fit result is shown in Fig. E^b). We find 15 600 ib 200 signal events. 

To determine the reconstruction efficiencies, two samples of signal Monte Carlo events for the 
decays D*~^ Dfir^ and D*~^ Dfj, each of which consists of 30 000 events, have been used. 
The Monte Carlo events have been analyzed using the same procedure as for real data. By cal- 
culating the ratio of reconstructed to generated event numbers, we find efficiencies of e{Df7:^) = 
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Am [GeV/c^] Am [GeV/c^] 

Figure 3: L>*° signals: (a) m{K-iT+ir°) - m(K"7r+); (b) m(K"7r+7) - m(K"7r+). The dots 

represent data points. The sohd hne shows the fitted function. The dashed hne indicates the 
portion of the fit associated with the background. 



0.041 ± 0.002 and e{D+-y) = 0.071 ± 0.002 for the two D*+ decay modes. The efficiency ratio is 
e(L>+7r°)/e(L»+7) = 0.58 ± 0.03. 

For completeness, and as a check of our efficiency calculations, we have also measured the ratio 
T{D*^ D^7r^)/T{D*^ —>■ D^^) using the same selection criteria for the vr*^ and photon candidates 
as in the reconstruction of D*'^ Dfir^ and D*'^ — > Df^y. To reconstruct the decay K~Tr~^, 
combinations of and vr^ candidates are required to successfully fit to a common vertex, and 
the scaled momentum of the resulting candidate must be 0.6 or greater. Fitting the sum of 
a double Gaussian function and a third-order polynomial to the resulting K~Tr~^ invariant mass 
distribution, we find (996.0 it 1.5) x 10^ signal events. K~tt~^ combinations are retained if their 
mass differs by less than 17MeV/c^ from the nominal mass jEj- 

The candidates are combined with all ir^ candidates; the resulting mass difference Am{D^7r^) = 
m{K~7r^7r^) — m(K~'ir'^) is shown in Fig. ISfa). A fit using a double Gaussian for the signal and 
the function shown in Eq. ^for the background yields 69 000 it 400 signal events. 

The candidates are then combined with all photon candidates producing the distribution 
of the mass difference /S.m{D^^) = m^K^ir^j) — m{K^iT^) shown in Fig. |21^b). In this case the 
signal peak of D*^ D^^ is close to a large bump due to the reflection of D*^ D^ir^ where one 
photon comes from a vr*^ decay. We therefore model the background using the function 

h{Am) = AT • 1^1 + exp ^"^ "J^^"""^ ^ ^ • (Am^ + aAm + b) . (3) 

(Note that this function is similar to Eq. ^ but differs in the sign of the exponential term.) 

The signal is modeled by a "Crystal Ball" function (Eq. [2)). The resulting fitted signal consists 
of 67900 ±700 events. 
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Table 1: Summary of systematic uncertainties on the branching ratio 
T{D*+ ^ Dt7T°)/T{D*+ ^ D+7). 



Signal Monte Carlo Statistics 5.0% 

Monte Carlo Simulation Uncertainty 3.6% 

Df Background 3.2% 

vr*^ Background 3.1% 

Momentum Dependence 6.8% 

Total 10.2% 



Table 2: Summary of systematic uncertainties on the branching ratio 



Signal Monte Carlo Statistics 


5.4% 


Monte Carlo Simulation Uncertainty 


3.8% 


Background 


0.1% 


Momentum Dependence 


2.8% 


Total 


7.2% 



We determine the efficiency ratio by the reconstruction of 30 000 signal Monte Carlo events for 
each of the channels D*° D^-k^ and L>*° D°j. Efficiencies of e(L'°7r°) = 0.037 ± 0.002 and 
e{D^-y) = 0.064 ± 0.002, and an efficiency ratio of e(L>°7r°)/e(L>°7) = 0.58 ± 0.03 are found. This 
is in perfect agreement with the value of e(L'+7r°)/e(L>+7) = 0.58 ± 0.03. A ratio of T{D*^ 
D'^n^)/T{D*^ ^ L>°7) = 1.740 ± 0.020 (statistical error only) is obtained, which agrees with the 
world average of T{D*° D°7r°)/r(D*° D'^-f) = 1.625 ± 0.200 [S]. 

Various sources of systematic uncertainties have been studied; they are summarized in Tables ^ 
and El 

To verify that the Monte Carlo events model the data correctly, we have studied r decays 
with one or two vr*^ mesons in the final state to obtain energy-dependent Monte Carlo efficiency 
correction functions for vr*^ mesons and photons. While no correction is necessary, the errors on the 
correction functions represent uncertainties of the Monte Carlo model and are used as systematic 
uncertainties. 

To test for uncertainties in the background shape of the mass difference distributions, we have 
considered left and right sidebands in the Df, D^, and vr" masses. The widths of the sidebands have 
been chosen such that the number of events in each sideband corresponds to the expected number 
of background events in the signal region. We have fitted the same functions used to determine the 
signal yields to the mass difference distributions of the sideband samples. Any discrepancy in the 
yield between data and Monte Carlo simulation is considered a systematic uncertainty. 

We have repeated the measurements of T{D*~^ D+7r°)/r(Z)*+ D+j) and r(D*° 
D^n^)/T{D*^ D^^) for candidates restricted to bins oi p* (encompassing the range ?> < p* < 
5GeV/c) by repeating the mass difference fit to the subsample of candidates within each bin. By 
fitting both a constant function and a first-order polynomial to each branching ratio as a function 
of p* , we have verified that the measured branching ratios are independent of p* . Nevertheless, 
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Table 3: Summary of preliminary results. The first errors are statistical, the second represent 
systematic uncertainties. 



T{D*+ - 


Dt7r')/T{Dr ■ 




0.0621 ± 0.0049 ± 0.0063 


B{D*+ - 






0.0585 ± 0.0043 ± 0.0056 


B{D*+ - 






0.9415 ± 0.0043 ± 0.0056 


T{D*^ - 






1.740 ±0.020 ±0.125 


B{D*^ - 






0.6351 ± 0.0027 ± 0.0166 


B{D*^ - 






0.3649 ± 0.0027 ± 0.0166 



we assume conservatively that the first-order polynomials arise from unknown momentum depen- 
dencies of the efficiencies which do not cancel in the branching ratios. Taking the momentum 
distributions of the D*'^ and D*^ mesons into account, we use the difference between the branching 
ratio represented by the constant function and the first-order polynomial for each distribution as a 
systematic uncertainty. 

The branching ratios are shown in Table El By assuming that the D^^ meson decays only to 
D+vrO and Z)+7, i.e. B{D*+ D+n^) + B{D*+ D+j) = 1, and that the meson decays 
only to L>°7r° and D^-f, i.e. B{D*^ L>°7r°) ± B{D*° D^j) = 1, it is possible to calculate the 
branching fractions, which are also listed in Table |31 

In summary, we have measured a preliminary value for the branching ratio T{D*~^ Dfir^)/ 
T{D*^ Df-f) = 0.0621 ± 0.0049(stat) ± 0.0063(syst), which is consistent with the previous 
measurement i2j, but has higher precision. We have also measured T{D*^ — > D^7r^)/T(D*^ — > 
L>°7) = 1.740 ± 0.020(stat) ± 0.125(syst), which is more precise than, and in agreement with, the 
world average J^. 

The branching fraction B{D*+ D+ti^) = 0.0585 ± 0.0043(stat) ± 0.0056(syst) is larger than 
the theoretical prediction, which is in the 1 — 3% range pp. However, it should be noted that the 
theory predicts a strong correlation between B{D*g^ Df-K^) and B{D*~^ — > D'^'y); for B{D*~^ — >■ 
< 1%, B{D*~^ Df-K^) is expected to be greatly enhanced. 
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